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Rank-based  Inference  for  Linear  Models 


Errors 


In  this  paper  robust,  rank-based  inference  procedures  are  considered 
for  general  linear  models  with  (possibly)  asymmetric  errors.  Approximating 
standard  errors  of  estimates  and  testing  hypotheses  about  the  model 
parameters  require  estimating  a  scaling  functional,  and  an  approach  is 
developed  which,  unlike  previous  work,  does  not  require  symmetry  of  the 
underlying  error  distribution  or  replicates  in  the  design  matrix.  Hence, 
important  asymmetric  models  such  as  arise  in  life  testing  can  now  be 
handled.  Further,  it  is  shown  that  the  asymptotic  properties  of  the 
inference  procedures  hold  with  simpler  conditions  on  the  design  matrix 
than  previously  required.  In  addition  an  estimate  of  the  intercept  is 
developed  without  requiring  the  assumption  of  a  symmetric  error  distribution 
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1.  Introduction  In  this  paper  robust,  rank-based  inference  procedures 
are  considered  for  general  linear  models  with  (possibly)  asymmetric 
errors.  Approximating  standard  errors  of  estimates  and  testing  hypotheses 
about  the  model  parameters  require  estimating  a  scaling  functional,  and  an 
approach  is  developed  which,  unlike  previous  work,  does  not  require 
symmetry  of  the  underlying  error  distribution  (McKean  and  Hettmansperger 
1976)  or  replicates  in  the  design  matrix  (Draper  1981).  Hence,  important 
;  asymmetric  models  such  as  arise  in  life  testing  can  now  be  handled. 

Further,  it  is  shown  that  the  asymptotic  properties  of  the  inference 
procedures  hold  with  simpler  conditions  on  the  design  matrix  than  pre¬ 
viously  required.  In  addition  an  estimate  of  the  intercept  is  developed 
without  requiring  the  assumption  of  a  synmetric  error  distribution. 

The  models  to  be  considered  and  the  basic  assumptions  are  now  given. 
The  vector  of  observations  ^  =  (Y1,  Y2,  Y3,  ...,  Yn)“  is  assumed  to 
satisfy  either 

(1.1)  X  *  +  *n  «  +  « 

or 

(1.2)  + 

where  ^  is  the  n  x  1  vector  of  ones,  ^  is  the  p  x  1  vector  of  unknown 
regression  coefficients,  a  -if  it  is  included-  is  the  unknown  intercept, 
and  =  (e^,  e2,  . .. ,  en) "  is  an  n  x  1  vector  of  independent,  identically 
distributed  randan  errors  with  continuous  cumulative  distribution 
function  F.  We  will  impose  the  following  assumptions  as  needed: 


(1.3) 


(1.4) 


F  has  density  f  having  finite  Fisher  information;  that  is, 
f  is  absolutely  continuous  and 


fcr 


(e)/f(e)r  f(e)de  <  •. 


Xn  is  a  known,  n  x  p  full-rank  matrix  with  ith  row 
The  centered  matrix  Cn  =  Xn  -  n"1  ^  ^  Xn  is  also  full 
rank. 


(1.5) 


(1.6) 


"  V'  On’  A  • 

a  positive  definite  matrix,  as  n  — >  00 . 

n"1  C;  C  — >  E , 
n  n 

a  positive  definite  matrix,  as  n  — >  00 . 


We  note  in  passing  that  (1,5)  implies  (1.6),  Where  possible  without  confusion, 
the  dependence  of  various  quantities  on  n  will  be  suppressed. 

The  inference  procedures  for  %  to  be  considered  are  based  on  a  measure  of 
dispersion  proposed  by  Jaeckel  (1972).  For  an  n  x  1  vector  define  the 


dispersion  as 
(1.7) 


DC#  s  a(i)  v(i)» 


where  v^  £  v^j  <  ...  <  are  the  ordered  elements  of  ^  and  a(l) ,. . .  ,a(n) 

are  a  set  of  scores  satisfying  seine  regularity  conditions.  Consideration  in 
this  paper  is  restricted  to  Wilcoxon  scores, 

(1.8)  a(i)  =  121/2  (i/Cn+1]  -  1/2). 

Thus  each  element  of  ^  is  assigned  a  weight  proportional  to  the  difference 
between  its  rank  among  the  n  elements  and  the  average  rank.  Procedures  based 
on  these  scores  generalize  the  Mann-Whitney-Wilcoxon  two-sample  procedure  and 
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inherit  its  asymptotic  efficiency. 

A 

Jaeckel  proposed  estimating  ^  by  minimizing  D(^  -  X@)  and  showed  the 

estimate  thus  obtained  is  equivalent  to  an  estiiiate  previously  suggested  by 

1/2 

Jureckova  (1971)  in  the  sense  that  n  times  the  difference  converges  to 
zero  in  probability.  Thus  the  two  estimates  have  the  same  limiting  distri¬ 
bution;  specifically,  under  conditions  to  be  weakened  in  this  paper  in  the 
case  of  Wilcoxon  scores,  n  converges  in  distribution  to  a 

multivariate  normal  random  variable  with  mean  Q  and  variance-covariance 
2  -1 

matrix  rl  ,  where 

(1.9)  t'1  s  121/2  r  f2(e)  de. 

J  —OO 

It  is  often  desirable  to  test  hypotheses  about  ^  of  the  form 

(1.10)  Hgj  =  Q  versus 

HA: 

where  H  is  a  full-rank  q  x  p  matrix  with  q  £  p.  A  consistent  test  for  (1.10) 
can  be  based  on  a  quadratic  form  in  the  full-model  estimate  of 

(1.11)  Q  =  T“2  frCH(C'C)"1  H']"1  H  £, 

A 

where  t  is  sane  consistent  estimate  of  t.  The  statistic  Q  has  an  asymptotic 
X  (q)  distribution  under  Hq,  If  q  <  p,  a  consistent  test  for  (1.10)  can 
also  be  obtained  by  fitting  both  the  full  model  and  the  reduced  model  induced 
by  Hq.  Letting  and  be  the  corresponding  estimates,  McKean  and 

Hettmansperger  (1976)  showed  that 

(1.12)  D*  =  2tDq  -  4^)  -  Dq  -  xifuu,):/^ 

2 

has  an  asymptotic  X  (q)  distribution  under  HQ  and  can  be  used  as  a  test 
statistic.  Both  test  statistics  require  a  consistent  estimate  of  the  scaling 
functional  t.  The  estimates  previously  proposed  necessitate  either  the 


assumption  of  a  synmetric  error  distribution  or  replicate  rows  in  X.  In 
Section  3  an  estimate  is  developed  without  either  of  these  assumptions. 

The  asymptotic  theory  for  Q,  and  D*  is  founded  on  the  asymptotic 
linearity  of  the  gradient  of  the  dispersion  D(^  -  x^)  treated  as  a 
function  of  jg.  Thus  the  technical  assumptions  adopted  by  Jureckova  (1971) 
in  her  proof  of  the  asymptotic  linearity  have  been  carried  over  by 
subsequent  authors,  In  Section  4  it  is  shown  that  this  linearity  property  - 
for  the  important  case  of  Wilcoxon  scores  -  can  be  obtained  without  seme  of 
the  complicated  assumptions  on  the  design  matrix  required  by  Jureckova; 
specifically,  her  assumptions  3a,  3b,  and  3c  are  eliminated.  The  results 
of  Kraft  and  van  Eeden  (.1972)  for  linearized  rank  statistics  should  also 
hold  under  less  complicated  assumptions  on  the  design. 

Since  DQ£  -  a^.  -  X£)  =  DQ£  X£)  ,  Jaeckel  ’  s  dispersion  function  pro¬ 
vides  no  information  concerning  the  intercept.  McKean  and  Hettmansperger 
(1978)  showed  that  a  can  be  estimated  by  applying  a  one-sample  signed-rank 
procedure  to  the  residuals  after  estimating  if  the  error  distribution 
is  synmetric ,  In  Section  5  an  estimate  is  proposed  which  does  not  require 

A 

symmetry  of  F  and  its  joint  asymptotic  distribution  with  £  is  stated. 

2,  A  Preliminary  Lemma  The  proofs  of  the  results  in  this  paper 
rely  on  a  lemma  which  is,  in  essence,  imbedded  in  the  proof  of  Theorem  3.1 
of  Jureckova  (.1969),  We  state  the  lemma  here  for  convenient  reference. 

In  the  sections  to  come,  we  are  concerned  with  the  asymptotic  behavior 

A 

of  some  random  variable  ^) ,  where  ^  is  a  consistent  estimate  of  a 

parameter  %  and  ^  is  a  random  vector.  For  simplicity  the  dependence  on 
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n  is  suppressed.  In  many  cases  it  is  easy  to  determine  the  behavior  of 

A 

and  if  it  can  be  shown  that  -  H(^r,^)  converges  to  zero  in 

A 

probability,  the  behavior  of  H(^,^)  is  also  determined. 

A  A  A 

For  example,  let  W  be  a  function  of  the  residuals  e^  =  ^  ^ 

after  fitting  a  linear  model.  Then  one  can  think  of  W  as  H((a  Q')',  y) , 
where  HC  (a  Jg ')  * ,  p  is  that  same  function  of  the  "residuals"  Y^  -  a  -  ^  Jg. 
Since  H((a  ,  %)  is  a  function  of  the  independent,  indentically  distri¬ 
buted  errors  e^  =  Y^  -  a  -  its  behavior  rtay  be  simple  to  determine. 

The  following  lenma  gives  sufficient  conditions  for  H(^,^)  -  H(^,^) 
to  converge  to  zero  in  probability. 

Lenma  2.1  Suppose  =  Q(g,^)  +  ^ (g-^f), where  Q(g,^)  is  monotone 

A«  A.  A»  A * 

in  each  of  the  components  of  g  ,  perhaps  nondecreasing  in  some, 

*  t 

nonincreasing  in  ethers  and  |  j£|  |/n  <  K  <  •  for  some  t  >  0.  Here  £  may 

depend  on  n;  |  j  •  1 1  is  the  supt-norra  throughout  this  paper.  If  H(g,^)  - 

% 

H(^,^)  converges  to  zero  in  probability  for  g  =  +  ^/n1"  with  ^  fixed, 

*>/ 

then  for  each  B,  0  <  B  <  ®, 

sup  |H(g,^)  -  H<^)| 

converges  to  zero  in  probability,  where  the  supremum  is  over 

(g:  n  I  |g  -  ;^|  |  <_  3h  Furthermore ,  if  n  is  bounded  in  probability, 

A>  Ai 

A 

then  H(^,^)  -  H(^,^)  converges  to  zero  in  probability. 

3,  Scaling  Functional  In  the  one-sample  setting,  several  authors 
(Bhattacharyya  and  Roussas  1969;  Schuster  1974;  Schweder  1975;  Ahmad  1976; 
Cheng  and  Serf  ling  1981)  have  considered  estimating  using  window  (kernel) 
density  estimates,  sometimes  as  a  particular  case  of  more  general  estimation 


problem.  However,  no  one  has  treated  the  problem  when  the  estimate  is 
computed  using  dependent  quantities ,  such  as  the  residuals  in  regression. 

Let  fn(e)  be  the  window  estinate  of  f(e)  given  by 

-1  n 

(3.1)  fn(e)  =  1  wCCe-e^J/h^) , 

where  e. ,  . . , ,  en  are  a  random  sample  from  a  distribution  having  density 
f,  w  is  a  density,  and  h^  Cn=l,  2,  3,  ,,.)  is  a  sequence  of  constants 
converging  to  zero.  Under  appropriate  regularity  conditions,  the 
estimate  for  6  =  t^/121^2,  given  by 

(3.2)  [  f  (e)dF  (e), 

J-oo  n  n 

where  Fn  is  the  usual  empirical  distribution  function,  is  strongly  consistent 
and  asymptotically  normal. 

The  estinate  (3.2)  can  be  written  as 

(3.3)  (n2h  )’'1  Z  Z  w(Ce.-e-]/h_) 

n  i=l  j=l  1  ]  * 

=  w(0)(rihn)“1  +  (n2!^)”1  ZZ _  w([ei-ej]/hn) . 

Since  the  nonrandan  contribution  of  the  i=j  terms  is  of  smaller  order  than 
the  randan  partial  due  to  the  i*j  terms,  for  the  asymptotic  theory  we  work 
only  with  the  latter.  Thus  we  consider 

(3.4)  i  =  [n2^]"1  IX  wtCe^-ej]/]-^) , 

A 

and  show  that  6  remains  consistent  for  6  when  the  e^  are  replaced  by  the 
dependent  residuals  e^  =  Y^-& -  ^^(or  e^  =  -  ^^)  and  hn  is  replaced  by 

A  A  A 

a  randan  h  ,  under  mild  conditions  on  w,  h  ,  fi,  and  h  .  In  applications  it 
is  necessary  to  use  the  data  to  determine  the  window  width,  in  order  to 

A 

obtain  good  performance  and  to  make  9  scale  equivariant.  Note  that  in 


rj-rp^Ti 


this  section  ^  is  not  necessarily  the  rank  estimate. 

Theorem  3,1  Suppose; 

(i)  w  is  a  square -integrable,  strongly  unimodal  density,  synroetric 
about  zero,  with  finite  second  moment; 

(ii)  hn  =  n”rh,  where  0  ^  r  <_  1/2  and  h  is  a  positive  constant; 

(iii)  n  (^-^)  is  bounded  in  probability;  and 

(iv)  n^(h*ph)  is  bounded  in  probability  for  sane  q  >  0 ,  and  hn  =  n“rh 
Then  under  assumptions  Cl. 3),  Cl. 4),  and  Cl. 6),  for  model  Cl. 1)  with 


ei  s  Yi  "  a  “  *Si£  (or  for  110(361  (1*2)  ei  s  Yi  “  $£§>> 

C3.5)  6  =  Cn2hnJ’"1  II  wCCe^-e^]/}^)  converges  in  probability  to  6. 


C3.6) 


Proof;  For  arbitrary  p  x  1  Jg  and  k  >  0  define 

TCb;k)  =  Cn2hn]-1rtw(Ki-^Jj-Y.^J/[n-rk]) 

=  CnV1 


where  =  n"*rk.  Then  TC,g;h)  is  the  estimate  of  6  based  on  the  independent, 
identically  distributed  e^  and  on  the  nonrandom  window  width  h^;  and 

^  A  ^  A 

0  =  Chn/hn3TC^;h)  is  the  estimate  based  on  the  residuals  e^  and  random  window 

A  A  A 

width  hn.  It  is  now  shown  that  T(^;h)  -  T(^;h>  converges  to  zero  in 


probability.  Since,  under  the  conditions  of  the  theorem,  T(£;  h)  is  easily 
shown  to  converge  to  6  in  probability  (Aubuchon(1982)  used  a  projection 
argument  to  show  that  n  (T(^;h)  -  0)  converges  in  distribution  to  a  normal 
random  variable  under  these  conditions)  and  hn/hn  converges  to  1  in 
probability,  this  implies  that  0  also  converges  to  0  in  probability. 

Although  TC£;k)  is  not  of  the  form  required  by  Lemma  2,1,  it  is 
possible  to  split  TQj>;k)  into  two  pieces,  each  of  which  is  of  the  necessary 
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Let 


form.  To  this  end,  we  define  the  following  monotone  functions. 

^3*7)  w1(z)  =  w(z)  ,  if  2  <_  0 

=  w(0)  ,  if  z  >  0;  and 

,  if  z  £  0 

w(z)-w(0)  ,  if  z  >  o. 


w2(z)  =  o 


Then  w(z)  =  w^z)  +  w2(z).  Also,  let  the  2p  x  1  vector 


(3.8) 


where  is  the  vector  of  positive  parts  of  ($£-i$j)  and  is 

the  vector  of  negative  parts,  so  that  )+  +  )~;  recall 

that  is  the  1th  rcw  of  the  design.  Letting  (far  2p  x  1 


(3.9) 


T*c£ Sk)  «Cn\]‘X  a  wde^j-JJy  >3/^), 


:  *13 

we  have  T*(£;k)  =  T(Jg,k)  when  £*=  (Jjf  %)•  Finally,  note  T*(^;k)  =  T*(^;k)  +  T*(^;k) 
where 
(3.10) 


T*<4;k)  =  Cn2^]'1  a  CV^i-«j-ei'j('«-(I))1/kn) 


-C-l)m0w(O) 

for  m  =  1,  2. 

The  following  lertma  establishes  conditions  under  which 
T*(^:k)  -  T* ( (^ "  ^")  ';h)  converges  to  zero  in  probability  for  m  =  1,  2.  The 
proof  of  the  lensna  involves  lengthy,  calculations  and  is  deferred  to  Appendix 
A.  The  linear  term  subtracted  in  T|  and  added  in  Tjj|  plays  a  crucial  role 
in  the  proof. 

Lenina  3 . 1  For  m  =  1 ,  2 ,  under  the  conditions  of  Theorem  3.1, 


laa 


“  n  -  where  4  =  C$'  fr  .  ^2  ^  k  ,  h  +  a/nq  f<jr 


Note  that 
(3.11) 


any  fixed  ^  and  a. 


U.ll)  Tm(^;k)  ”  Tm((|)’h) 

=  tn2hnrli£/.-»m<Cei-ej-«iV«-(|»^n> 

-(-uVy1:!  e»co)  n{^-(p). 

The  first  summation  is  monotone  in  the  components  of  while 


2,  ,-l 


II Cn*hn]  x  EE.  ^ij  1 1  /h"*"7  ^  £.  2(n*  hn)-^  m^x  n”^E  Cl+x^^],  vrfiich  is  bounded. 

L  A 

Since  n  ($-^)  is  assumed  to  be  bounded  in  probability,  Lenina  2.1  (with  t=l/2) 
implies  that  c 


(3.12) 


T*(U);k)  -  T*(C|);h) 


converges  to  zero  in  probability  for  m  =  1,  2  and  k  =  h  +  a/r.0-.  Furthermore  , 

A 

g  Q  A 

T*((-');k)  -  T|5EC('S);h)  is  monotone  in  k.  Since  nH(h-h)  is  assumed  to  be 
m  a  ’  m  p 

* 

bounded  in  probability,  it  follows  from  Lenina  2.1  that 


(3.13) 


S  „  Q 

T*(U);h)  -  T*(0§);h) 
m  Q  7  mb 

£  * 


converges  to  zero  in  probability  for  m  =  1,  2, 


Recalling  that 
(3.14) 


T(#h)  -  T(#h) 

=  T*(A;h)  -  T*C($);h> 

S  * 

1  -  ft 

=  T  *((.);h)  -V((«);h) 

1  3  ~  1  % 

*  &  ,  a 

+  T,ft(Q;h)  -T„*((«)-;h), 
1  S  « 


this  yields  the  desired  result. 


% 


’■  %  ‘.V.  A* 


(3.15) 


xantple  For  applications,  we  consider  a  modified  form  of  (3.3): 

b*  =  (nh)"1  +  [n(n-l)h  T1  II  w([e.-e, ]/h  ) , 

n  i/j  1  3  n 

where  h  is  a  constant,  and  we  take  h  =  n”^^h.  The  authors  (1984,  Section  4) 

n 

show  that,  in  the  iid  case,  the  bias  is  then  of  order  n*1.  If  the  underlying 
error  distribution  has  density  f(y,6)  =  <5~^f^(6~^y) ,  then 


(3.16) 


h  =  C2'1 


[f{(y)r  dy 


u^  w(u)  du]-^^ 


will  make  the  first  -  order  terms  in  the  bias  vanish. 

If  we  take  f^  to  be  the  normal  density  with  interquartile  range 

A 

equal  to  1,  6  the  sample  interquartile  range  (defining  h) ,  and  w  the  uniform 
density  on  (-1/2,  1/2),  then  by  Theorem  3.1, 


(3.17)  0*  =  (4.11n  3)_1  + 

is  a  consistent  estimator  of  6. 


( 4 . lln17  2 (n-1 ) 6 ) -1  II  w 


4.  Linearity  of  the  Gradient  of  the  Dispersion  In  this  section,  it 
is  shown  that,  for  Wilcoxon  scores,  the  asymptotic  linearity  property  of 
the  gradient  of  the  dispersion,  first  proved  by  Jureckova  (1971),  holds 
under  simpler  conditions  on  the  design  matrix.  Thus  the  work  of  subsequent 


authors  in  developing  the  asymptotic  properties  of  Q,  and  D*,  which 

relied  on  the  linearity,  also  is  valid  under  simpler  assumptions. 

The  dispersion  of  ^  -  X  may  be  written  as 

n 


(4.1) 


I 

i=l 


a(R(Yi^))(Yi^) 

Z  CRCfi^/dt+l)  -  1/2](Y.-^), 


=  12 


1/2 


where  ^  is  the  ith  row  of  the  centered  design  matrix  C  (1.5)  and  R(Y^-jg^Jg) 
is  the  rank  of  Y^-jgjjg  among  Y^-*g^Jg,  . . . ,  Then  the  gradient  of 

D(^-X  Jg),  taken  as  a  function  of  Jg,  exists  except  on  a  set  of  Lebesque 
measure  zero  in  tiP.  When  it  exists,  the  negative  of  the  k  element  of 
this  gradient  is 

(4.2)  n1/2Sk(^)  =  121/2  cikCR(Yi^)/(n*l)  -  1/2] 

=  121/2  (n+l.-1  RW.^Jj). 

Since  the  ranks  in  (4,2)  are  translation  invariant,  S^Cjg)  can  be  further 
rewritten: 

(4.3)  S^)  =  121/2Cn1/2(n41)]-1  Cik  R(eiM£r(jj-|)) , 

noting  that  et  =  Y±  -  a  -  3^.  (or  e£  =  Yi  -  $$) . 

Theorem  4.1  For  model  (1.1)  ar  model  (1.2),  under  assumptions  (1.3), 

(1.4) ,  and  (1.6),  and  far  S^(Jg)  defined  by  (4.2), 

(4.4)  supl^Jj)  -  S^)  +  T-1  JJ(k)n1/2()j-^)  | 

converges  to  zero  in  probability  as  n  -*■  ®,  where  the  supremum  is  over 
n1/2  |  |  <_  B  <  oo  f  and  is  the  kth  column  of  E,  the  limit  of  n^C'C. 

Proof  Slutsky's  Theorem  implies  that  (n+lf^Qg^ ),  where  is  the 

k  column  of  C,  may  be  substituted  for  since  (n+1)  C'C  -*Ias  n  **■  00 . 
Define 

(4.5)  Tk(J6>  =  Sfcfly  4  T-1  (Q§(k))-  n1/2  (Jj-^)/(n41). 

Again  we  cannot  directly  apply  Lenina  2.1  to  T^(Jg)  but  must  split  the 


quantity  into  two  pieces.  As  before  (3.8)  let 


(4.6) 

'«  •  (?J-) 

\  Cj£i  3  / 

,  /  (Si-s/ 

Letting 

\  (?i-?j)' 

(4.7) 

Tkl(^)  =  121/2[n1/2 

-1  n 

(n+1)]  1  Z  C. 

n 


-  T“1[n(n+l)]_1n1/2(  Z  cj.  I  D..)'(i-  (&)  and 

i=l  lk  *  « 


Tk2(#  = 


121/2[n1/2(n+l)]-1 


n 


n 


C'iD  £  nej-ei<^i(^(|):} 


-  x’1Cn(n+l)]‘1n1/2(  ?  CT.  S  D-. ) '  (£-(& ) 

i=l  **  j?l  3  ^ 


for  any  2p  x  1  we  have 

(4.8)  Tk(Jj)  =  Tj^q)  4  T^q)  when  £  *  (Jj'  Jj')' 


The  following  lenma  establishes  that  Tj^C^)  -T^C  ^  ')  ')  converges  to 
zero  in  probability  for  m  =  1,  2  when  ^  =  (£"  ,§'')"  +  ^/n  .  The  proof 

of  the  lemma  involves  lengthy  calculations  and  is  left  to  Appendix  B. 


Lemna  4.1  For  m  =  1,  2,  under  the  conditions  of  Theorem  4.1, 

ecV0-ti«  *  o 

as  n  -*>  «  where  ^  Q')  '  +  j^/n^2  for  any  fixed  d.  Now  the  first 

portions  of  T^($)  and  T^^)  are  monotone  in  the  components  of  Futher, 

1 1  Cn3/2(n+l)]“1n1/2  Z  C*  Z  |  and  ||[n3/2(n +l)r1n1/2  Z  C“  Z 

i=l  1K  j=l  -11  i=l  j=l  J 

are  bounded.  Thus  Lemma  2.1  implies  that 

(4.9)  supit^-t^ct  r>'>! 


converges  to  zero  in  probability,  where  the  supremum  is  over 


n1/2  114-  Cg'  j6')'II  _<  B>.  Now  A  =  n1/2  |  -  <£'  £'>'|  |  <  B 

and  =  (£'  Jg')}  is  a  subset  of  this  set.  Thus  (4.9)  still  converges  to 
zero  in  probability  when  the  supremum  is  taken  over  A.  This,  along  with 
(4.8),  yields  the  desired  result,  since  the  supremum  of  a  sum  is  less  than 
or  equal  to  the  sum  of  the  supreme . 

A 

5.  Intercept  A  simple  estimate  a  of  the  intercept  for  model  (1.1) 

A  A 

is  given  by  the  median  of  ,  Yn~^.  Using  a  proof  similar  to 

that  of  Theorem  4.2  (c)  in  McKean  and  Hettmansperger  (1978),  it  is  possible 
to  show  that  n  7((a-a)  (£-£) ') '  converges  in  distribution  to  a  multi¬ 
variate  normal  random  varible  with  mean  Q  and  variance-covariance  matrix 

'(4  f2(0))_1  +  tVx  -t2ux  E"1 

V=  -t2  r\  *2ri  J, 

where  Z  is  the  limit  of  n-1  C'C  and  j£x  is  the  limit  of  the  p  x  1  vector 
of  column  means  of  X,  The  conditions  needed  are  (1.3),  (1.4),  and  (1.5); 
See  Aubuchon  (1982)  far  the  details. 


Before  proceeding  with  the  proof,  some  preliminary  consequences  of 
the  assumptions  of  Theorem  3.1  are  stated  without  proof. 

1/2 

Lamia  A.  1  For  any  ^  of  length  2p,  ^/n  -*•  0  uniformly  in  i,  j 

as  n  -*■  00  . 

Lama  A.  2  f  is  bounded. 


Lama  A.  3  If  G  is  the  cumulative  distribution  function  of  e^  -  e2> 
where  e^  and  e^  are  independent  and  have  c.d.f.  F,  then 

(i)  the  density  of  G, 

g(z)  =  [  f(e+z)f(e)de, 

J  —CD 

is  bounded  and  absolutely  continuous; 


(ii)  the  derivative  of  g(z). 


g'(z)  -  j  f"(e+z)f(e)de, 

J  -00 

is  bounded  and  absolutely  continuous;  and 

(iii)  g"(z) ,  the  derivative  of  g’(z),  is  bounded. 

Lama  A.  4  For  any  ^  of  length  2p,  n-2£E  (J^d)2  is  bounded. 

Lonna  3.1  For  m  =  1,  2, 

ECT*C^;k)  -T*((|);h)]2  -  0, 

where  £  =  (.%'  +  ^/n1/,2and  k  =h+a/n^,  as  n  -*■  ®,  for  any  fixed  ^  and  a 

Proof  Consider  the  case  of  m  =  1.  Referring  back  to  (3.10), 

ECT|C,$;k)  -  T*((£'  h)]2 


=  n“V’2  ll  ll  ECCw(e.-e.-J3'.4/n1/2]/k JKe^e.+^/n172) 
n  i/j  s^t  1  3  13  n  3  3 

-  w([ei-e.]/hn)I(ei<ej} 


(A.  1) 


% 


1 


i 


! 


V.* 

;< 

■ 


♦  w(0)(I{ei>ej+^/n1/2}-I{ei>ej  }+6^j4/n1/2)  ] 

xCw([es-et-^/nl/2]/kn)I{esiet+J2ij^/nl/2} 


-w(Ces-€t3/htl)I{eslet} 


+wCO)  CIies>et+J^/n1/2  }-I{es>et  l+e^/n172 )  ]] 
Partition  the  terms  in  this  quadruple  sum  into  three  groups: 


Group 

Description 

Count 

Gl 

Two  matching  pairs  of  subscripts 

2n(n-l) 

G2 

One  Hatching  pair 

4n(n-l)(n-2) 

G3 

Mo  Hatching  pair 

n(n-l)(n-2)(n-3) 

n2(n-l)2 

TOTAL 

and  deal  with  each  group  separately. 

Consider  the  terns  in  group  Gl;  they  are  uniformly  bounded  (using 
Lemma  A.l),  so  their,  sum  is  2n(n-l)n~4hn”20(l)  =  (nhn)“20(l)  0 

as  n-*-00,  where  0(1)  is  bounded. 

Among  the  terms  in  group  G2  consider,  for  example ,  the  sum  of  those 
with  j  =  t: 


(A.3)  |n-\-2  £  My  x  My  f(8j)  d6j  | , 

j=t 

1  "-V2  |2  Jl  '"ij!  x  |Msjl  f(ej>  de: 

3*t 

<  n"4h‘2  I  sup | Mj • |  sup|M  • | , 

“  n  G2  3  J 

j* 

where 


15 


■ 


(A.  4) 


fe>y/n1,t  1/2 

Mij  =  J_.  w([ei'e:'Sij^/n  ]/kn  >  f<ei)dei 

t  ei 

-  w([e.-e.]/h  )f(e. )de. 

i  ]  n  1  i 

+  w(0)|"  Cl{ei>e.+^j4/n1/2}-I{ei>ej}+6^/n1/2]f(ei)de. 

and  the  suprema  are  over  all  i,  j ,  s ,  and  .  Now  if  sup | |  -*•  o  and 
sup|M  .  |  -*■  0  as  n  00 ,  then  the  sum  of  terms  in  group  G2  with  j  =  t 

SJ 

3  2  1 

coverges  to  zero,  since  there  are  0  (n  )  terms  in  the  sum  and  (nh  )  is 
bounded. 

But  letting  u  =  [e^-e j  2  ]  /  kn  in  the  first  integral  and 

u  =  [e.-e. ]/h  in  the  second, 
a.  ]  n 

(A-5)  lMijl  s  lkn  w(u)f(ej+knu+^^/n1/2)du 

f° 

"hn  J  w^JfCej+h^uJdu 

4w(0)CF(e.)-F(ej+QCj^/n1/2)+e^rj^/n1/2]| 

iCKn+hp)  sup^f^2 

+wC0)[6  +  sup(ffl/,|  %'j^/n17  2 1  -  0 

uniformly  in  i  and  j. 

Similar  arguments  hold  for  the  other  cases:  i  =  s,  i  =  t,  and 
j  =  s. 

New  turn  to  the  sum  of  terms  in  group  G3. 


(A. 6) 


n’4h  "2  Z  E[( •)(•)] 
n  G3 

=  n"V“2  Z  EC-]  E[-], 


since  the  multiplicands  are  independent  when  no  subscripts  natch. 
Consider  one  of  the  expectations: 


(A.  7) 


E[w([ei-e-^rjd/n1/2]/kn)I{ei-ej  <  J^/n372} 

^wCCe^^jl/l^Jlfe^-ej  <  Q> 
twlIlKHe^,  >JJy4/n1/2)-I{ei^,  >  0} 

♦  9)JyiJ/n1/2)] 


rP'-d/n1'2  -  /0 

=  j  1D  w([z-J2^/n1/2]/kn)gCz)dz 


r° 

w(z/hn)g(z)dz 

* 


+w(0)[GC0)-GC]gCj^/n^/2)  +  6^.^/n172], 

recalling  that  G  is  the  cdf  of  e^-e^  with  density  g,  when  i*j .  Letting 
1/9 

u  =  [z-^j^/n  ]/kn  in  the  first  integral  and  u  =  z/h^  in  the  second , 
this  expectation  is 

(A. 8)  kn  f  w(u)g(knu  +  £T^/n1/2)du 

J  *00 


fO 

-hn  wCuJgCh^  u)du 
•  *00 

-wdJXJ&jJ/n1'2)2  G"(ii5)/2, 


i  •  V»  Vv.  -  ./v' vS  >/.’  vi/.v ' 


wrsew 


where  |C  —  |  <  I ,  recalling  that  6  =  G'(0).  Making  a  Taylor 

expansion  of  g( • )  about  zero  in  each  integral,  further  reexpress  the 
expectation : 

<A'9)  knf  w<u^g^J+(V+^ij^/nl/2)8'(0)+(knu+J3ij^/nl/2)2g"(i;ij)/2]du 

“hnj  0  nug  ^  CO )  ■*n2u2g"  ( <. )  /  2  3du-w(  0 )  (^j^/n^  2 )  2G"  (£  ^ ^ )  /  2 , 


where  <  |knu+^j^/n1/2 1  and  j 4 j  S  |knu| .  Now  g'(0)  =  0  and  g"(‘) 

r°  r° 

is  bounded  (Lenraa  A.3).  Further,  w(u)du  =  1/2,  |  ju|w(u)du«»,  and 
•0 


_oo 


r»: 


w(u)du«».  Thus  we  can  write 


(A.  10) 


n"2h  “1  £Z  I  EC*  31 
n  i*j 


’  "'V1  «  t| Vhnl*(0)/2+lknl3cltl)tkn2l«yJ3/"1/2l0<1) 


i^j 


+  lkn  ^/nl/2)2|0(1>V0(1)+(^/nl/Z)2°(1):i  -  0 


1/2,2, 


ay 


as  n-v«,  since  l^-hj/l^  *  °»  1^1 3/hn  -  °»  kn2/hn  -  °>  I$ij3/nl/2 1 

_2  2  *3 

uniformly  in  i,  j,  n  ZZ.  ^ij$)  ^  bounded  (Lama  A.4)  and  (nh^)  -►  0. 


But 


I  n”4h”2  l  ECOEC'31  <  Cn'V1  ZZ  |E[-3|32  -  0  as 
n  G3  -  n  i*j 


n-+>x>.  Therefore  the  sum  of  terms  in  group  G3  converges  to  zero.  A  similar 
argument  holds  for  T*(£;k)  -  Tjj|((|);h),  and  the  proof  of  the  lerrma  is 
complete. 
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Appendix  B:  Proof  of  Lenina  4.1 

Before  proceeding  with  the  proof,  a  preliminary  lemma  which  follows 
fran  the  conditions  of  Theorem  4.1  is  stated  without  proof. 

-1/2 

Lenina  B.l  n  max  |C..  |  *►  0  as  n  -►  ®,  for  k  =  1,  2,  ...,  p. 

l<i<n  ** 

Lenina  4.1  Far  ms  1,  2,  under  the  conditions  of  Theorem  4.1, 

-  vcr  ° 

1/2 

as  n  +  ®,  where  6  =  (Q'  Q')  '  +  j£/n  for  any  fixed 

Proof  Far  simpler  notation,  the  factor  (n+1)  in  T^C^)  is  replaced  by 
n;  see  (4.7). 


Proof:  Calculating  the  expectation  far  m  =  1, 


(B.l) 


ect^)  -  t^' 4'nj2 


n 


n 


s  12n-3  ECjc^  jFiCl(ej-ei  <  ^  ^/n1/2> 
-Kej-e^  <  0}  -  e^/a1/2 


X  s=l  t=lCI{*t‘*s  1  ’J/nl/2> 


1/2- 


-I{et-es  <  0}  -  6^s  $/n^]] 


s  12n"3  ZZ  ZZ  ctv  cL  EC  (*)(•)  3 , 
i*j  sft  **  SK 


since  =  0.  The  terms  in  this  quadruple  summation  can  be  partitioned 


into  two  groups: 

(B.2) 


Group 

<31 

G2 

TOTAL 


Description 

Some  matching  subscripts 
No  matches 


Count 

4n(n-l)(n-2)  +2n(n-l) 
n(n-l)(n-2)(n-3) 

7 - 


7 


n  (n-l)‘ 
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(B.3) 


Consider  the  sum  of  terms  in  the  first  group: 

12n‘3  £ 

=  12n'3  £  4  4  MCI(ej-ei<I2j'i4/n1/2> 

-  1{ej^i  £  °>3  *  ^1{et^s  £^ts  ^1/2}-I{et-es  <  0}]] 

*  12n‘3  G\  4  4  a-2%#«v«s  i  ^"1/2> 

-  Het-es  <  0}-  6)Jjs4/n1/2]/n1/2 

-  ®2 

*  Qjl  *  q2. 

Now 

(B.4)  |QJ  <.  12n-3  £  <4  <4  EC| -|  x|-|3 

GX 

<  12T3  £  44  CE| -|  xEl-H172, 

G2 

by  the  Cauchy  -  Schwarz  inequality,  noting  that  1 1{  •}  -  I{  •} j  is  equal 
to  its  own  square.  Furthermore, 

<B.5)  ElKe.^  <  ^/nl/2}-  I  (e^  <  0}| 

*  |G%jJ/n1/2)  -  G(0)| 

where  |£j<J  £  l^ji^nl/2|*  ^  G'  (Lenina  A. 3),  and 

I/O 

|jg4j^/n  |  converges  to  zero  uniformly  in  (i,j)  as  n  -*•  (Lenina  A.l). 
Further, 
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(B.6) 


"~3  a  ^  ** 


•1  ft  •  ^  o  .  . 

<  n*1  Z  (c  )z  +  n~*  E  E  c+  c+ 

i=l  **  i=l  3=1  ^  sk 

■i  n  ^  -I  n  - 

<  n"1  E  +  [n  1  E  |c.J]  , 


i=l 


i=l 


•ik1 


n 


£  ^  .*  <4  *  Cl  +  H*1  *  CL]2 

i=l  **  i=l  ** 

which  is  bounded  since  n-1  C'C  converges.  Thus  the  right-hand  side 
of  CB.4)  converges  to  zero  as  n  and  Q1  does  also. 

Now  consider: 


(B.7)  Q2  =  12n-3  £  4  <4  [-^[GCJJ^/n1'2) 

G1  J 

-  G<0>  -  ^34/n1/J]/n1/2  -  62<Jjr.(j)(ftt'34)/n] 

--  22"'3  2  4  4  [-24<3<43/nl/2)2  G"  etsV(2n1/2) 

G1 

where  l5tgl  <  Ij^gfJ/n222! ,  recalling  that  G^CO)  -  j  f2(e)de  -  6  . 

Since  G"  is  bounded  (Lenina  A. 3),  n  E  c^.  c^  is  bounded  (B.6), 

G1 

>  1/9 

and  |  -*■  0  uniformly  in  (i,j),  we  have  Q2  *►  0  also.  Referring 

back  to  (B,3),  the  contribution  of  terms  in  group  1  converges  to  zero 


as  n  +  «. 

Finally,  turn  to  the  second  group.  These  terms  must  go  to  zero 
faster  than  n”1,  since  there  are  0(n4)  of  them  and  the  divisor  in  front 

O 

is  only  n  .  Since  no  subscripts  match,  the  multiplicands  are  independent 
and  the  expectation  factors: 
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Thus  by  showing  that  the  sum  over  terms  in  each  of  the  two  groups 
converges  to  zero,  we  have  shown  that  the  expectation  on  the  left-hand 
side  of  (B.l)  converges  to  zero.  This  establishes  the  lenma  for  m  =  1 
The  proof  for  m  =  2  is  analogous. 


'•.VJ'T  V  I"*  ^  ’.1  V*.  I.\  .V  t’A'.-.V.'  -.■ 
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parameters  require  estimating  a  scaling  functional,  and  an  approach  is 
developed  which,  unlike  previous  work,  does  not  require  symmetry  of  the 
underlying  error  distribution  or  replicates  in  the  design  matrix.  Hence, 
important  asymmetric  models  such  as  arise  in  life  testing  can  now  be 
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handled.  Further,  it  is  shown  that  the  asymptotic  properties  of  the 
inference  procedures  hold  with  simpler  conditions  on  the  design  matrix 
than  previously  required.  In  addition  an  estimate  of  the  intercept  is 
developed  without  requiring  the  assumption  of  a  symmetric  error  distribution. 
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